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Abstract 

The low-temperature thermal conductivity in polycrystalline graphene is theoretically studied. 
The contributions from three branches of acoustic phonons are calculated by taking into account 
scattering on sample borders, point defects and grain boundaries. Phonon scattering due to sample 
borders and grain boundaries is shown to result in a T a -behaviour in the thermal conductivity 
where a varies between 1 and 2. This behaviour is found to be more pronounced for nanosized 
grain boundaries. 
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I. INTRODUCTION 



The thermal conductivity of both single-layer and few-layer graphene is of current ex- 
perimental and theoretical interest (see a recent review [l[ an d the references therein). It 
was shown that the heat conduction in graphene has a specific behaviour due to the unique 
nature of two-dimensional phonons. Main attention has been attracted to studies of the heat 
transfer in graphene near room temperature because of possible applications in electronics 
and photonics. 



Our interest is stimulated 
typically polycrystalline {2- 



In this paper, we focus on the thermal conductivity in graphene at low temperatures. 

yy recent observations that the large-scale graphene films are 
5] and consist of many single-crystalline grains separated by 
grain boundaries (GB) [6|-[8|]. The grain sizes are dependent on growth conditions, ranging 
from hundreds of nanometres to tens of micrometres for slight changes in growth condi- 
tions 6j. Recently, the fabrication of nano crystalline graphene made by electron-radiation 
induced cross-linking of aromatic self-assembled monolayers and their subsequent anneal- 



ing was reported 



9]. Nanosized grain boundary loops originating from paired five- and 



seven-membered ring disclinations were observed and analysed in 10]. 



The GB-induced phonon scattering has a marked impact on the thermal conductivity, k, 
at low temperatures. In bulk materials, the well-pronounced crossover from k ~ T 2 to k ~ 
T 3 has been proved with the crossover temperature strongly depending on the GB size 



121 ] . Therefore, one can expect a similar effect in two-dimensional polycrystalline materials 



like graphene. The paper is organized as follows. In Section 2 we overview the general 
formalism with three main acoustic phonon branches taken into account: longitudinal (LA), 
transverse (TA) and flexural (ZA). The effective temperatures characterizing the phonon 
mean free paths for main scattering mechanisms are introduced. The results of numerical 
calculations of the thermal conductivity are presented in section 3. Conclusion is devoted 
to the discussion of the results obtained. 
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II. GENERAL FORMALISM 



According to [13j, the phonon thermal conductivity of single-layer graphene can be 
written as 

K = . . T2h £ / (nuJ s (q)) 2 v s (q)l s (q, T) x 
;T 2 h f ,f f „ Jo 



4irk B T 2 h e ff 



e hui a (q)/{k B T) 



where summation is performed over the phonon polarization branches with the dispersion 
relations u s (q) (q is the wavevector, qp> corresponds to the edge of the Brillouin zone, and 
s is the branch index), h e ff is the effective graphene layer thickness, v s (q) = du s (q)/dq is 
the phonon group velocity, l s (q,T) is the phonon mean free path, and kp is the Boltzmann 
constant. The phonon mean free path, l s , is limited by three major independent scattering 
mechanisms: sample borders (rough boundary (RB)), point defect (PD), and grain bound- 
ary. Within the relaxation-time approximation, it reads 

l/l s = I/Iq + 1/lpo + l/l GB , (2) 

where Iq, Ipo and Iqb comes from phonon scattering by RB, PD, and GB, respectively. The 



phonon-RB scattering is given by 



15| 



/-i _ I LlP ( 3 ) 



with p being the specularity parameter defined as a probability of specular scattering at the 
sample borders and d the graphene layer size. The phonon-PD scattering can be written as 

7-1 _ 5 oP q(u) 2 U) 

where Sq is the cross-section area per one atom of the lattice and V is the measure of the 
scattering strength. It was suggested in that grain boundaries can be naturally described 
by disclinations. In graphene, the alternating pentagon-heptagon structure along high-angle 
GBs was revealed in [8]. This finding allows us to model the grain boundary by biaxial 
wedge disclination dipole (BWDD) of length C which consists of 5-7 pairs (see Fig. [Q. In 



this case, the phonon-GB scattering can be expressed as [11] 



l GB = 2D 2 (uC) 2 n t qg(qC), (5) 



FIG. 1. The fragment of grain boundary in polycrystalline graphene. The size C of this wall 
exceeds eight lattice translations. Notice that the whole wall corresponds to a single BWDD. 



where v is the Frank index, rii is the areal density of BWDDs, 

Q{z) = J 2 (z) + Jl(z) - J (z)J 1 (z)/z, 
Ji(z) is the Bessel functions of ith kind, and 

D = 7T 7 (l-2<7)/(l-<7), (6) 

where 7 is the Gruneisen constant of a given phonon branch and a is the Poisson constant 



sec 
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for detail). 

Generally, there are six phonon branches. However, at low temperatures only three 
branches are of importance: acoustic longitudinal (LA) and transversal (TA) phonons with 
the ordinary dispersion relation u s = v s q (s = LA, TA) and out-of-plane (flexural) acoustic 
phonon branch (ZA) with uza = q 2 /2m, where m = 2 J P2d/K is an effective parameter fl^ . 



P2D arid K are graphene 2D mass density and bending stifness, respectively. Notice that at 
sufficiently low temperatures the umklapp processes of phonon-phonon scattering are mainly 
frozen and can be neglected. 

Let us introduce the dimensionless parameter x, x — fuoj {k B T). For LA and TA phonons 



eq.(pQ) takes the form 

C q T 2 r°D/T _ x 3 e x dx 



Ktf Jo x ^ T) (^Ty> (7) 

where 9 D is the Debye temperature, C s = k B l /(4:iih 2 v s ), and 

x;\x,t) = 1 + (x^-) 3 + x^-g(x^) (8) 

V J-pdS j-gb j-q 
is the total mean free path normalized to Iq (A s = l a /lo)- As a matter of convenience we 
have introduced here three characteristic temperatures: 

for PD scattering, 

To = 7 r^, (10) 

and Tg_b = hv s /(2D 2 i> 2 C 2 kBlo n i) for GB scattering. In order to clarify the role of the GB 
scattering let us temporary disregard the second term in eq. (jSJ) which is responsible for PD 
scattering. In this case, I, ss / (l + 2T /(7rT Gi j) + O(T 3 /(T Gi? T 2 ))- 1 for xT » T (the short 
wavelength limit). As is shown below, typically To > Tgb- Therefore, namely To defines the 
region of marked influence of phonon-GB scattering and can be considered as a threshold 
temperature. It should be also mentioned that To characterizes the temperature when the 
wavelength of an incident phonon becomes comparable with the size of grain boundary. 

To take into account the imperfect packing of grains one can introduce the average 
distance between grain boundaries, a, so that the areal density of grains is estimated as 
rii = (a + £)~ 2 and 

_ hv s (l + a/£) 2 
TgB ~ 2D 2 u 2 k B l ' (H) 
where (1 + a/C)~ 2 is a packing coefficient. Notice that while T GB depends on a/C, T is a 

function of £ only. 

The heat conductivity due to ZA phonons is written as 

C ZA T 3 / 2 fOzA/T x 5 / 2 e x dx 

k za — r / ^za{x,T)— — — , (12) 

h eff Jo [e x - iy 

5/2, 



where C ZA = y/2mk 5 J 2 l / (47ih 3/2 ) and 



X z 1 Ax,T) = l+ 1 2 ZA (x,T)((x^-) 3 / 2 + 

-LPD 



+ f£w- (13) 



The characteristic temperature of PD scattering takes the form 



rro = ±(«£f» (14) 

m \S Tl J k B 
while for GB scattering one obtains T = h/(2C 2 mkB) and 

rpt Hi + a/cy 

GB ~ 8mk B [D'^H }^ 

where D' = tt(1 — 2a)/(l — a). Notice that the Gruneisen constant for ZA phonons 7(g) 
is assumed to depend on the wavevector q in contrast to LA/TA cases. In accordance with 

n 

[1 91 ] one can use an approximate expression 7za(<?) = — 1 — ^{q/qd — l) 2 or, equivalently, 



1za(x,T) = -1 - S0(Jx^- - l) 2 . (16) 

V Oza 



III. RESULTS 



The total phonon thermal conductivity of graphene is calculated by using eqs. ([7]) and 
()12p . The major initial parameters used in the model are d = 5/im, p = 0.9, Vla = 21.3km/ s, 
vta = 13.Qkm/s (from [13]), 7la = 1-8, 7ta = 0.75 (from |l|), and jza(^,T) determined by 
eq. f fT6l) (which approximates the dependence of Gruneisen constant for ZA branch on the 
wavevector), as well as m = 1.0 ■ 10 6 s/m 2 are taken from One should note, that the 

Gruneisen parameter for the ZA phonon branch is variating within the Brilluin zone by the 
order of several magnitudes, in contrast with LA and TA cases. To take into account the 
influence of ZA phonons, we use eq. (Tl6|) for the Gruneisen parameter instead of constant 
value. We use the upper limit for the specular graphene thickness parameter h e ff = 0.69nm 



from 



3- 



One should also note that the differences in the recent theoretical and experimental 



models give rise to the discrepancy in the determination of certain graphene constants (see 



the discussion in |l| for detail). However, this discrepancy will not play crucial role for the 
effects discussed below. 

Using the constants defined above, one obtains Cla — 8.40 ■ 10~ n W/K 3 , Cta = 1-31 • 
10- W W/K 3 , C Z A = 5.93 -10- 9 W/K 5 / 2 . The characteristic temperatures are given in tables H] 
and HIl We restrict our consideration to the temperatures below 100 K. To estimate the 
parameter Tpd we fit our model in the case of graphene single crystal according to the 
results of [3]. We found that T PD exceeds 100i^ for the leading phonon modes, and thus PD 
scattering turns out to be of little consequence within the temperature interval considered. 




1 2 5 10 20 50 100 1 2 5 10 20 50 100 

T, K T, K 



FIG. 2. The total thermal conductivity of polycrystalline graphene with high-angle grain bound- 
aries for C = 2A6nm (left) and C = 7.38nm (right) at different packing coefficients: (a) a/C = 0, 
(b) a/C = 2, (c) a/C = 9, (d) a/C = 30, and (e) a/C = 100 (see table HI|) . Dashed curve cor- 
responds to a case of graphene single crystal. Dot-dashed straight lines shows T 1 ' 5 and T. The 
approximate crossover temperatures are indicated with vertical lines. 



FigfS] shows the thermal conductivity of graphene as a function of temperature for C = 
2.46 nm and £ = 7.38 nm at different packing. For comparison, the dot-dashed lines shows 
the functions T L5 and T, while the dashed curve corresponds to a case of graphene single 
crystal. The thermal conductivity of the graphene single crystal follows a square law in 
a wide temperature region excepting the very low temperatures where kza has a marked 
impact and high temperatures (in the range of 80 — 100 K) where PD scattering has a slight 
influence. 

For polycrystalline graphene, one can clearly see the characteristic points where crossover 
occurs. The crossover temperature is estimated as T* w T /2, where T is governed by the 
dominant phonon mode, (TA). Above T* the thermal conductivity behaves like T 2 while 
below T* the behaviour is acutely sensitive to Tqb (see eq. ffTTj) ). In turn, Tqb depends 
distinctly on the packing coefficient and GB's angle. Notice that high-angle GBs with 



v ~ 0.083 were revealed 



in 



For this case, the typical values of Tqb at different packing 
coefficients are shown in table [III At low Tqb one obtains k ~ T while growing Tqb results 
in k ~ T 15 and, finally, in k ~ T 2 . One can conclude that the lower is Tqb, the more 
pronounced influence of grain boundaries takes place. 



0.05 
0.02 
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FIG. 3. Reduced thermal conductivity, k/T 2 , versus temperature T for a graphene single crystal 
(dashed) and polycrystalline graphene (solid) at fixed packing coefficient a/C = 9 and different 
grain sizes: (a) C = 2A6nm, (b) C = 7.38nm, and (c) C = 22.14nm. Dot-dashed straight lines 
shows T 1 ' 5 and T. 

TABLE I. The threshold temperature at different grain sizes for three phonon branches. 

C, nm 2.46 7.38 22.14 
T ,K (LA) 66.0 22.0 7.33 
T ,K (TA) 42.1 14.0 4.68 
T ,K (ZA) 0.6 0.066 0.0074 



Fig. |3] shows reduced thermal conductivity, k/T 2 , versus temperature at fixed packing 
coefficient. As is clearly seen, both the thermal conductivity and the crossover temperature 
decreases with increasing GB size. In case of graphene single crystal, one can mention once 
again a slight deviation from a T 2 -behaviour due to the influence of ZA phonons. Point 
defects are responsible for deviations at high temperatures. 
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TABLE II. The grain boundary temperature at different effective distances between grains for the 
case of high-angle grain boundaries (u = 1/12). 



a/C 


0.0 


2.0 


9.0 


30.0 100.0 


T GB , mK (LA) 


1.5 


13.5 


149.7 


1438.9 15273.6 


T GB , mK (TA) 


5.5 


49.5 


550.6 


5291.8 56172.5 


T' GB ,mK (ZA) 


1.3 • 10" 


' 7 1.06 • 10~ 


" 5 0.0013 


0.12 13.63 



IV. CONCLUSION 

In this paper, the influence of polycrystalline structure on the low-temperature heat trans- 
fer in suspended graphene has been theoretically studied. Three acoustic phonon branches 
were taken into account and three types of scattering mechanisms were discussed: sample 
boundary scattering, point defects and phonon scattering on the grain boundaries. The 
rough-boundary scattering is found to play an important role over the whole considered 
temperature range (0-100 K). On the contrary, the point defect scattering manifests it- 
self only at the upper limit of the temperature interval. The influence of polycrystalline 
structure is described by two temperature parameters Tqb & n d T which, in turn, define 
the crossover temperature T*. Above T*, the heat conductivity behaves like in graphene 
single crystal, n ~ T 2 , with an adjusted rough-boundary scattering parameter. Below the 
crossover point we found that k ~ T a where a values range from 1 to 2 with increasing 
Tqb- 111 experiment, such behaviour of the thermal conductivity could be a manifestation 
of the grain-boundary scattering mechanism. Tqb decreases with an increase of the packing 
coefficient, the misorientation angle and/or a sample size while T depends on the grain 
boundary length only. 

In our study, the contribution of TA phonons is found to be dominant. The contribution 
of LA phonon branch is suppressed due to the higher (in comparison to TA) sound velocity 
while ZA phonons are of little importance due to the large Gruneisen constant. To be more 
precise, in our consideration ZA phonons manifests themselves only at very low temperatures 
and the grain-boundary scattering additionally suppresses their influence. 



It is interesting to note that in some recent experiments with suspended single-layer 20 



2 if ] and few- layer 



22] graphene the temperature dependence of k was found to follow a 



power law with an exponent of 1.4± 0.1. This behaviour can be attributed to a dominant 

9 



contribution from the ZA phonon modes. However, the measured thermal conductivity in 
those samples are much lower compared to the theoretically expected values. As possible 
alternative explanations one considers scattering of phonons in the bilayer graphene by 
a residual polymeric layer, a high concentration of defects due to processing [23|, or the 
polycrystalline structure of graphene with small and misoriented grains [24J. In this article 
we have shown that the phonon-GB scattering will objectively result in T a behaviour of 
thermal conductivity with 1 < a < 2. 

This work has been supported by the Russian Foundation for Basic Research under grant 
No. 12-02-01081. 
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